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ABSTRACT 

Large scale, weakly collimated outflows are very common in galaxies with large infrared luminosities. 

In complex systems in particular, where intense star formation (SF) coexists with an active galactic nu¬ 
cleus (AGN), it is not clear yet from observations whether the SF, the AGN, or both are driving these 
outflows. Accreting supermassive black holes (SMBHs) are expected to influence their host galaxies 
through kinetic and radiative feedback processes, but in a Seyfert galaxy where the energy emitted in 
the nuclear region is comparable to that of the body of the galaxy, it is possible that stellar activity 
is also playing a key role in these processes. In order to achieve a better understanding of the mech¬ 
anisms driving the gas evolution specially at the nuclear regions of these galaxies, we have performed 
high-resolution three-dimensional hydrodynamical simulations with radiative cooling considering the 
feedback from both star formation regions including supernova (type I and II) explosions and an AGN 
jet emerging from the central region of the active spiral galaxy. We computed the gas mass lost by 
the system, separating the role of each of these injection energy sources on the galaxy evolution and 
found that at scales within one kiloparsec an outflow can be generally established considering intense 
nuclear star formation only. The jet alone is unable to drive a massive gas outflow, although it can 
sporadically drag and accelerate clumps of the underlying outflow to very high velocities. 

Keywords: galaxies: Seyfert - galaxies: ISM - galaxies: outflow - galaxies: AGN -ISM: kinematics 
and dynamics - ISM: evolution 


1. INTRODUCTION 

Gas outflows extending for few kiloparsec scales (kpc) 
or larger are frequently observed in local galaxies with 
high star formation rates (SFRs) and in active galactic 
nuclei (AGN). A large fraction of the outflows is detected 
in high-redshift, active star forming sources, and for this 
reason it is common to believe that outf lows represent 
a sta ge of the galactic evolution (see, e.g.. iShaplev et ahl 
l2QQ3f ). Such outflows affect the multiph ase distributi on of 
the interstellar medium (ISM; see, e.g.. lMcKe3lI995[ ) and 
are responsible for the dynamical and chemical evolu¬ 
tion of the galaxy, determining the enrichment of the in- 
tergalactic rnedinm (IGM; see e.g. iHeckman et al .111 9901: 
iMelioli et al.ll2QQ9L I2QI3LI2QI5L and references therein) as 
well as the deficit of baryons seen in many galaxies. 

Outflows may also quench star formation by heating 
up the cold gas co mponent and ejecting it from the host 
galaxy (see, e.g.. iBinnevI I2QQ3 ). High energy outflows 
have been proposed as the main cause of the drop in AGN 
luminosity and the end of the galactic gas fuelling onto 
the central super massive black hole (SMBH), limiting 
the active phase of the nucleus to a time of about 10^ 
yr. Outflows are also particularly important for starburst 
(SB) galaxies whose intrinsically large SFR provides the 
perfect conditions for gas outflow to develop. 

In complex systems where an intense star forma¬ 
tion (SF) activity coexists with an AGN, it is unclear 
whether the SF, the AGN, or both are driving the out¬ 
flows. It is known that accreting SMBHs influence their 
host galaxies through kinetic/radiative feedback pro- 
cess es, but the details of these processes remain unclear 

(e.g. iStorchi-Bergmann et al.l[2QQ9L l2QIQl : iMorganti et al.l 
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I2QI3I ). In particular, in Seyfert galaxies, where the nu¬ 
clear source emits an amount of energy comparable to 
that of the host galaxy, it would be possible that the 
outflows are driven predominantly by the stellar activity, 
which therefore would have a key role in the evolution of 
the system. 

There are several studies that assume a ther¬ 
mal AG N feedback in t he v i cinity of the black 
hole (e.g JdT Matteo et al.l I2QQ5I: iSpr ingel et al.l I2QQ5I: 
iJohansson et al.ll2QQ9[ ). Observations indicate th at there 
may be two distinct modes of AGN feedback (see lFabi^ 
l2QI2l . for a review). The kinetic-mode (also denominated 
radio mode), thought to occur at low Eddington ratios 
(L/I/Edd < 0-1), is expected to operate through heat in¬ 
jection into the surrounding gas by a relativistic jet. In 
the quasar-mode, likely linked to high Eddington ratios 
(I//I/Edd > O.I), the AGN power is believed to couple 
directly to the interstellar medium of the host galaxy via 
radiation pressure or accretion disc winds, resulting in 
fast outflows (e.g. lAalto et al.ll2Q12l: iGicone et ^l2Ql4 

iGenzel et al.l l2Ql4 iTombesi et al.l l2QI5f ) . The general 
conclusions of these studies is that AGN feedback can un¬ 
bind a large fraction of the baryons provided that a few 
percent of the bolometric luminosity couples thermally 
to the ISM. This applies specially in the framew ork of 
quasar or radiative-mode AGN feedback (see. e.g..lKing 
2003 : lOstriker et al.ll20I0l: iFancher-Gignere fc Qnataert 

2niliZiihovas fc Kiridl2m2n 

IGaspari et alJ ()2011bi yi and lGaspari et, alJ (l2fll2D per- 
formed three-dimensional simulations of the evolution of 
the gas in clusters, groups of galaxies and elliptical galax¬ 
ies taking into account the AGN feedback and found that 
a mechanical feedback with these lower efficiencies would 
be enough to both constrain the growth of the black hole 
and balance the radiative cooling resulting in gas den- 
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sity and temperature distributions compatible with the 
obse rvation and with the results of si milar studies (see, 
e.g.. iFalceta-Goucalves et alJl2QlQbl lal: iCiotti et alJl2QlQl 

iNovak et al.ll2Qlffr 

However, these studies, in general, have focussed on the 
large (100 kpc — 1 Mpc) scale evolution of the gas, where 
most of the energy of the SMBH jet is directly transferred 
to the intra-group or intra-cluster medium carrying out 
a large amount of gas and probably preventing cooling 
flows. Almost no work has addressed the effects of the 
AGN on the host galaxy, where the highly collimated jet 
has possibly very weak interactions with the surrounding 
medium. Of course, there is a theoretical interest in the 
ability of AGNs to affect the gas mass i n the nuclear 
region and its star formation history (e.g.. iGroton et aTI 
l2QQ6l : IS^awinski et al.ll2QQ^ . but the inclusion of effects 
like the interplay between star formation and the AGN 
activity, the merger of galaxies and secular processes, in 
general, is complex, and the parameter space is not yet 
fully mapped, mainly at small (1 kp c ) scal es. 

In a recent work, I Wagner et al.l (j2Q12f ) investigated 
with help of HD simulations the detailed physics of the 
feedback mechanism between a SMBH jet and a two 
phase fractal ISM in the kpc-scale core of a galaxy, but 
their results refer only to the early stages of the jet evo¬ 
lution, when it propagates through the disk for the first 
time. In this case the dense embedded clouds of the ISM 
are accelerated by the ram pressure of the high velocity 
flow through the porous channels of the warm phase gas, 
but probably this regime cannot be maintained for much 
time after the jet digs a hole in the disk and beyond it and 
therefore, this cannot be taken as a stationary solution. 

Generally, feedback of the ISM implies not only en¬ 
ergy and m omentum injection by sn pernovae and stellar 
winds (e.g. iMcKee fc Ostrikerl[l97^ . b ut also radiation 
from massive stars and from the SMBH (jThompson et al.l 
l2QQ5l . and references therein). The coupling between 
photons and the interstellar gas is expected to be very 
weak if only Thompso n cross-section is implied (e.g. 
iGiotti fc Ostrikerll2QQlf ). but this coupling is enhanced 
in the presence of dust. Dust grains are able to scat¬ 
ter and absorb UV radiation re-emitting it into infrared. 
Thus since the gas is coupled to the grains hydrodynami- 
cally, radiation pressure on the dust may be able to push 
the gas against its self-gravity providing an additional 
driving to gas outflow. 

When momentum from stellar radiation pressure and 
radiation pressure on larger scales via the light that 
escapes from star- f ormin g regions is considered (see 
iHopkins et al.l[2QTTL |2Q12[ ). a more realistic multi-phase 
ISM develops, maintaining a reasonable fraction of the 
ISM at densities where the thermal heating from super¬ 
novae has a larger effect. The potential importance of ra¬ 
diative momentum fe edback has been explored for some 
time now (see, e.g., iHaehneltl 119951: IGiotti fc QstrikeJ 
l2QQ7f ). and its potential role i n the context o f AGN and 
SBs was first po i nted out by iMurrav et al.l (j2QQ5f ) and 
iThompson et al.l (j2QQ5f ). It might be considered in prin¬ 
ciple at least in some situations to study the evolution of 
the gas of a galaxy hosting an AGN. 

In this study we will consider kinetic (or radio) mode 
AGN feedback and ne glect the radiation and dust drive n 
mechanisms (see also (MartiniI2QQ5I: iSharma et al.ll2Qllf ). 
The photo-ionizing radiation will be considered only to 


justify the minimum temperature of the diffuse gas, 
Tmin,^ ^ 10^ K. In fact, we assume that all the me¬ 
chanical energy and momentum are effectively deposited 
into the ISM, so that the mechanical feedback will domi¬ 
nate over the radiative feedback. This is a reasonable 
assumption for Seyfert galaxies for which the column 
densiti es are about ^ 10 ^^ cm^ (or less), as pointed 
out by I Novak et al.l (|2Qllf ). Moreover, although dusty 
structures, in the form of spiral and filaments at hun¬ 
dreds of parsecs scales ar e often observed in early- type 
active galaxies (see, e.g., iSimoes Lopes et ^ 120071 ) in¬ 
dicating that a reservoir of dust is a necessary condi¬ 
tion for the nuclear activity, in Seyfert galaxies the dust 
mass is lower than that observed in ultraluminous IR 
galaxies (ULIRGs), and a rapid destruction in supernova¬ 
generated shock waves may make the typical dust lifetime 
shorter than the time over which the galaxy is active 
(IMnKepinflSa rT ones fc Niit,Ml2(ril]) . 

In order to achieve a better understanding of the pro¬ 
cesses driving the nuclear gas evolution in a Seyfert 
galaxy, in this study we perform fully three-dimensional 
(3D) hydro dynamical simulations at high (1.9 pc) reso¬ 
lution taking into account the effects of radiative cooling 
and considering the feedback from both star formation 
regions with type I and H SNe and from a collimated 
SMBH jet propagating from the central region of an ac¬ 
tive spiral galaxy. We compute the gas mass lost by 
the system and separate the role of each energy injection 
source on the galaxy evolution. 

The paper is organized as follows. In Section 2 we sum¬ 
marize some recent observational results about the gas 
dynamics in the inner kiloparsec scale of Seyfert galaxies. 
In Section 3 we describe the main physical processes in 
the nuclear region of an active galaxy and present solu¬ 
tions for a steady state system. In Section 4 we outline 
the main characteristics of our model, while in Section 5 
we present the main results of the simulations and finally, 
in Section 6, we summarize our conclusions. 

2. OBSERVATIONS 

Seyfert galaxies are usually spiral galaxies whose nu¬ 
clei are exce ptionally brigh t. The original definition 
of the class ((Sevfertl Il943l ) was primarily morpholog¬ 
ical, i.e., these are galaxies with high surface bright¬ 
ness nuclei, but today the definition has evolved so that 
Seyfert galaxies are now identified spectroscopically by 
the presence of strong, high-io nization emission lines. 
iKhachikian fc Weedmanl (|1974l ) were the first to real¬ 
ize that exist two distinct subclasses of Seyfert galaxies. 
Type I and Type H, which are distinguished by the pres¬ 
ence or not of broad bases of permitted emission lines. 
The spectra of Type I Seyfert galaxies show broad lines 
that include both allowed lines, like HI, Hel or Hell and 
narrower forbidden lines, like OHI {broad line regions^ 
BLR), while the spectra of Type H Seyfert galaxies show 
only (permitted and forbidden) narrow lines {narrow 
line regions^ NLR). In some cases the spectra show 
both broad and narrow permitted lines, and these ob¬ 
jects are then classified as an intermediate type between 
Type I and Type H, such as Type 1.5 Seyfert. NLRs are 
characteristic of a low-density ionized gas with widths 
corresponding to velocities of several hundred kilome¬ 
tres per second, while BLR have widths of up to 10"^ 
km s“^, and the absence of broad forbidden-line emis- 
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sion indicates that the broad-line gas is of high density, 
so the non—electric—dipole transitions are collisionally 
suppressed. BLRs are associated to the nuclear region 
< 1 pc) where the relativistic SMBH jet develops, 
while NLRs are observed between 10 pc and 1 kpc, and 
are associated to less dense clouds that may be optically 
thick and/or optically thin. 

The physical conditions of the NLRs have been ob¬ 
served accurately, and in the last years a large num¬ 
ber of studies have shown that photoionization by 
the central continuum s ource plays a dominant role 
in exciting the gas fe.g. iKraemer fc Crensha^ l2QQQal: 
iKraemer et al.l I2QQQI: iBennert et al.l l2QQ6blla[ ). although 
shocks may play an important r o le in localized re¬ 
gions (IKraemer fc Crensha"^ l2QQQbl: iDopita et al.l |2QQ3: 
iGroves et al.l 1200*^ Generally, the NLR is the largest 
observable feature which may be affected by the AGN 
radiation and its dynamical forces. Therefore, its inves¬ 
tigation is crucial to understand the AGN structure and 
evolution, as well as the interaction between the nuclear 
engine and the circumnuclear ISM and the stellar ac¬ 
tivity of the host galaxies. Despite a large number of 
ground-based long-slit spectroscopic studies, up to now 
a consensus about the kinematics of the NLR and its 
velocity distribution has not been reached yet. Several 
dynamical models h ave been proposed such as infall, ro¬ 
tation, outflow (e.g I VeillenxI 119911: iFraanelli et al.ll2QQQL 
l2QQ3l:[Fdscher et al.ll2QlQl:lMhller-Sanchez et al.ll2Qll[ and 

references therein), and radial acceleration by radio 
j ets or tangential e x pansion of the ga s around the jets 
(|Winge et al.l 119971: I Axon et~aLl Il998f ). However, ob¬ 
servations with the Space Telescope Imaging Spectro¬ 
graph indicated that the main component of the mo¬ 
tion in the NLR seems to be radial and more recently, 
further observational constraints on the kinematic of 
the NLR in Seyfert galaxie s have been established (e.g. 
iMiiller-Sanchez et al.l [2QTTI ). Though biconical models 
can explain very well the radial velocity profiles of the 
outflows, there are still several local variations in the gas 
kinematics that remain unexplained. 

Also recently, observations of Seyfert galaxies based 
on integral-field spectrosco py have found evidences of 
radial outflows in the N L R (iGarcfa-Lorenzo et al]l2QQll: 
Storchi-Bergmann et al.l I2Q1QI: iMhller- Sanchez et al.l 
20nl ). In Tese studies, it were measured the kinematics 
of the ionized and molecular gas surrounding the nucleus 
of the galaxy and it was found that the ionized-gas kine¬ 
matics may be the result of three velocity components, 
namely, an extended emission at a systemic velocity 
observed in a circular region around the nucleus, an 
outowing component along the bi-cone characterized by 
an angle of ^ 60°, and another component due to the 
interaction of the radio jet with the galactic disk. The 
data indicate that the origin of the extended emission is 
gas from the galactic plane that could be ionized by the 
AGN and suggest that the NLR clouds may be accel¬ 
erated very close to the nucleus 10 pc), after which 
the flow moves at essentially constant velocity. The 
observed mass outflow rates vary between 1 and 10 Mq 
yr^, being ^ 10^ to 10^ times larger than the expected 
accretion rates for an AGN, while the kinetic powers 
associated to the outflow are between 10^ and 10^ times 
smaller than the typical AGN b olometric luminosities 
(|Storchi-Bergmann et al.l 120101: iMhller-Sanchez et al.l 


Imll). Regarding the third kinematic component, that 
is, the one due to the interaction of the radio jet with 
the ambient gas moving at a systemic velocity, there are 
clear indications that the jet is launched close to the 
plane of the galaxy, pushing the gas it encounters on 
its way. The observed high-velocity gas in the nuclear 
regions of some galaxies as well as the large ratios 
between Mout and M^cc indicate that the outflow cannot 
be attributed to a single AGN wind episode. It might be 
rather due to a process in which the interstellar matter 
around the nuclear region is pushed aw ay by a previous 
AGN wind (|Muller-Sdnchez et al.ll2QlTI) . 

Molecular hydrogen emission arises in extended regions 
along the axis of the galaxy stellar bar, avoiding the re¬ 
gion of the bi-cone. The H 2 velocities are close to sys¬ 
temic with a small rotation component, supporting an 
origin in the galactic plane. On the contrary, the max¬ 
imum velocity of the outflowing component ranges be¬ 
tween 120 km s~^ and 2000 km s“^ and is reached at a 
typical distance of ^ 180 pc from the AGN. 

3. STEADY STATE SOLUTIONS 

As described, the physical processes occurring in the 
nuclear region of a Seyfert galaxy are quite complex. SNI 
and SNII explode at a given rate injecting energy, mo¬ 
mentum and metals into the ISM. The supernova rem¬ 
nants (SNRs) generated by the SN explosions expand 
through a turbulent multiphase medium in which clouds 
and filaments are embedded. At the same time, the gas 
heated by the SNe and by ionizing photons emitted by 
massive stars is also cooled via radiative losses that de¬ 
pend on the square of the gas density and on the metal 
abundance affected by the stellar activity. At the center 
of the galaxy, the presence of the SMBH also contributes 
to make more complex the evolution of the system. In 
fact, imaging studies have revealed that structures such 
as small-scale disks or nuclear bars and associated spiral 
arms are freque ntly observed in the inner kiloparsec of 
active galaxies ([Erwin fc Soar kelll9991: iPogge fc Martini! 
I2QQ2I: iLaine et al.ll2QQ3[) . transporting gas into the inner 
few hundred parsecs. In the same central region, small- 
scale jets emerge from the base of the accreting disk 
of the black hole, transfering a large amount of energy 
into a very collimated region, producing strong shocks 
and bright emission knots typical of collisionally excited 
plasma. As stressed in §1, in this contest it is very hard to 
understand if the observed gas outflow is a consequence 
of the SMBH feedback or if SNe-driven hot superbubbles 
may fill out the core of the galaxy developing a wind rich 
in filamentary structures and turbulences, or both effects 
are prevailing. 

3.1. Energy balance 

A first analytical approach study of the evolution of 
this sort of system may be done by comparing the energy 
injected by the main sources (SNI, SNII, Jet) with the en¬ 
ergy lost by radiative cooling and the energy required to 
drive a gas outflow. The rate of SNI explosions is about 
0.01(Mb„;ge,io) yr“^ (jPain et al.lll996[) . where Mb„;ge,io is 
the stellar mass of the bulge in units of 10^^ Mq, while 
the rate of SNII ex plosions correspon ds to about 0.01 
(SFRs'bt) fsee lMelioli et al.ll2009[) . where SFR^^,! 
is the star formation rate of a given SB region in units of 
1 Mq yr“^. Therefore, assuming that each SN releases 
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10^1 the injected luminosity associated to the SN 
events is: 


Lin = 3 X 10^^ [SFRsb,i + Mbhj] erg s ^ (1) 

where Mbrj is the mass of the central SMBH in units 
of 10^ Mq and where we have considered the observed 
relation between the bulge and the SMBH mass. 

The total rate of radiative cooling in an optically thin 
plasma within a cylindrical region of radius Rsb ,300 (in 
units of 300 pc) and thickness h 5'^^200 (in units of 200 
pc) is: 0 

Lcooi = 10^^ (A(T)_23 Rsb 300 hsB ,200 nfo) erg 

( 2 ) 

where nio is the gas density in units of 10 cm ^ at the 
height z = 0, and A(T )_23 is the cooling function in 
units of 10“^^ erg s“^ cm^, corresponding to a temper- 
atnre of about 1 0^ K. Given the Kennicutt-Schmidt law 
(|Kennicnttll 1 9981) that provides an excellent parametriza- 
tion of the global SFR over a density range extending 
from the most gas-poor spiral disks to the cores of the 
most luminous SB galaxies, after some algebraic opera¬ 
tions we obtain 

SFRkpc ~ 10“^ (nio hsB,20o)^'^ M©yr“^kpc“^ (3) 

where SFR/^pc = 3.5 SFR^^/R^qq, and substituting Eq. 

[3] in Eq. [2] we may directly compare the energy injected 
with the energy lost as a function only of the SER (per 
kpc^) of a given galaxy. Therefore, the rate of energy 
radiated away may be rewritten as: 

ieooJ = 2.7x10^3 (a(T)_23 RIb,soo *si ,200 SFR^^f) erg s 

( 4 ) 

and the trend of the energy injected and lost (per second) 
by the system is shown in Eig. [TJ 

Regarding the jet associated to the SMBH, an amount 
of energy Ebh ^ O.IMbhc^ could in principle be in¬ 
jected in the surrounding gas. This means that about 
bxlO^^MBHj erg s“^ may be released over a time of 10^ 
yr, a value significantly larger than the binding energy 
of the interstellar medium in a typical spiral or ellipti¬ 
cal galaxy, or even in a group of galaxies. However, in 
a Seyfert galaxy the luminosity associated to the central 
active jet is lower, between 10^^ and 10^^ erg s“^, i.e., of 
the same order of the luminosity injected by the SNe (Eq. 

[1]). It is also comparable to the energy carried out per 
second by the gas mass outflow, that is Eout ^ 

^ As remarked in the introduction we are neglecting dust driven 
mechanisms and radiation pressure because of the low dust column 
densities observed in the nuclear regions of the Seys and the ex¬ 
pected high frequency of SN shocks waves that will help to destroy 
dust structures. At the same time, the shocks as well as photoion¬ 
ization and non-thermal heating processes by the AGN will ensure 
a high gas temperature (we will find below in the simulations val¬ 
ues as large as 10® K; see. e.g., Figures 10 and 17). Thus, the 
assumption of an optically thin gas in near ionization equilibrium 
at the pc-scales around the AGN is reasonable. In this case, the 
thermal radiative cooling function given in Eq. 2 should be a good 
approximation even in the presence of photoionization and non- 
thermal heating by the AGN. This will actually help to keep the 
surrounding interstellar gas at temperatures larger than or equal 
to 10-* K. 



Log SFR (Mg yr-‘) 


Figure 1. Rate of energy injected and lost in the nuclear region 
of a Seyfert galaxy considering an active region with cylindrical 
geometry characterized by a radius R=300 pc and a thickness h= 
200 pc as a function of the SFR. Solid (black) line: energy injected 
by SNI and SNII; long-dashed (magenta) line: energy injected by 
the SNI only; dashed (blue) line: energy lost by radiative cooling 
assuming a fixed value for the cooling function A(T) = 10“23 erg 
s“^ cm“^; dot-dashed (red) line: energy lost by radiative cooling 
assuming a fixed value for the cooling function A(T) = 10“22 erg 
s“^ cm“^. The SFR is given in Mq yr“^ (logio scale) while the 
_epergy rates are given in units of 10"^^ erg s“^ (logio scale). 

Moreover, unlike the events described above, near the 
central black hole, the jet is strongly collimated and its 
interaction with the ISM of the host galaxy is limited to 
a small fraction of the total gas mass of the system, of the 
order of 10 “^(Rjet, 20 /^ 300)^5 where Rjet ,20 is the width 
of the jet within the galaxy disk in units of 20 pc. There¬ 
fore, we can expect that, except for the very early stages 
of evolution, the SMBH jet will influence only weakly the 
increase of the mass outflow, so that we can neglect its 
contribution in Eq. [H 

Erom Eig. [U we note that for a cooling function with 
a value of ^ 10“^^ erg s“^ cm^, which is a lower limit for 
a gas at T=10^ K, the rate of energy injected is larger 
than the rate of energy radiated away only for a SER 
smaller than ^0.1 Mq yr“^ kpc“^, and for a value of 
^ 10“22 g-i which is an upper limit for a gas 

at T=10^ K, the rate of energy injected is larger than 
the lost energy rate only for SER smaller than 10“^ Mq 
yr“^. In fact, according with the Kennicutt-Schmidt law, 
the SER increases slower than the square of the density, 
and therefore there is a critical gas density value above 
which an outflow cannot develop because of the radiative 
cooling, as indicated in Eig. [TJ 

3.2. Hot gas and filaments 

The result above is mainly due to the idealized sce¬ 
nario adopted with a homogeneous distribution of the gas 
which misrepresents the expected and observed mnlti- 
phase ISM (see, e.g., iMcKee fc Qstrik^ll977l : iDettmarl 
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llQQQtlHni et al.|[2QT^ in the disk of the galaxies. Indeed, 
when the stellar winds and the SN energy are injected 
into the SB region, two competing processes occur: on 
one hand the winds and the SNRs sweep the diffuse gas, 
thus lowering its mean density; on the other hand the 
interactions between the shells of the several SNRs de¬ 
termine the formation of dense and cold structures, gen¬ 
erating turbulences and a network of filaments which may 
coex ist with the high-temperature, low-density gas (see, 
e.g.. iMelioli fc de Gonveia Dal Pinol I2QQ4 iMelioli et al'l 
l2QQ5f ). Therefore, the importance of the radiative cool¬ 
ing depends on which of these two processes dominate. 
Besides, the radiative cooling is almost negligible in the 
hot bubbles generated inside the SNRs. 

The average distance crossed by a SNR before interact¬ 
ing with another one depends on its expansion velocity 
and thus on the density of the ISM, and also on the rate 
of SN explosions per volume. Assuming that a SNR ex¬ 
pands for a maximum ti me ^ yr 

(see equations 6 and 9 in IMelioli et al.l[2QQ6f ). during this 
period the number of SN explosions is: 

Ktaii - 5 X 10" mo Rloohioo (5) 

and the average distance between two SNe is: 

\sN ^ 13 pc (6) 

Therefore, when the average distance between two SNe 
is smaller than the stall radius of the SNRs, these may 
interact developing a low density ambient at high tem¬ 
perature coexisting with filamentary structures. Figure [2] 
illustrates the average distance (solid lines) and the stall 
radius (dotted lines) as a function of the density and 
temperature of the ISM. We note that if the Kennicutt- 
Schmidt law is valid, when the temperature is 10^ K or 
lower, SNRs always interact for ambient densities higher 
than ^0.1 cm“^. When the temperature is about 10^ 
K, each SN explosion may be considered as a single event 
and there is no SNRs interaction mainly because of the 
high thermal pressure of the ISM. 

If there is interaction, the successive generations of 
SNRs expand in a medium more and more rarefied, thus 
becoming adiabatic and eventually transferring most of 
their energy to the ISM; in this case the heating ef¬ 
ficiency (HE) of the SNe, that is, the fraction of the 
SN explosion energy that remains effectively stored in 
the ISM gas and is not radiated away, will be close to 
unity. On the other hand, in a situation where fila¬ 
ments, clouds and clumps are destroyed on short time- 
scales raising again very rapidly the ambient density, HE 
will be very small (for a detailed study o n the HE see, 
e.g.. IMelioli fc de Gonveia Dal Pinoll2QQ4D . Therefore, in 
a multiphase ambient medium characterized by diffuse 
low density gas and denser filaments, if we assume an es¬ 
cape velocity of the order of ^ y/2GMsys/hz^ the mass 
loss rate of the gas in the nuclear region will be approx¬ 
imately given by: 

.V HExE _ 

M ~ - (7) 

^ GMsys ^ ^ GMsys ^ ^ ^ 

where E and jC are the injected and radiated power, re¬ 
spectively. Considering an average total mass of the sys- 
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Figure 2. Average distance, Xsn, between two SNe and SNR stall 
radius, RstalU versus density n, for three different temperatures of 
the ISM, T=100 K (red lines), T^IO^ (black lines) and T=10® K 
(blue lines). Solid line: \sN\ dashed lines: Rgtall- The distances 
are given in pc (logio scale). 

tern of ^ 10^^ Mq and a height hz = I kpc, the outflow¬ 
ing bipolar gas mass rate, as a function of the global SER 
for different values of HE, is given in Eig. [3l The most 
striking result is that, for values of SER smaller than 
or approximately equal to 1 Mq yr“^, the SB activity 
drives an outflow smaller than 1 Mq yr“^, regardless of 
the SER. Eor SERs larger than 1 Mq yr“^, different sce¬ 
narios are possible. Depending on the efficiency of the 
radiative cooling, a strong outflow may develop or not. If 
HE is high, between 0.5 and 1, a mass loss rate between 
25% and 50% of the SER is possible, but for smaller HE 
values either an outflow or an inflow can occur, or even 
both of them concomitantly. Of course, the higher the 
SER, the density surface and/or the jet interaction with 
the ISM, more complex is the description of the evolu¬ 
tion of the whole system using a simple approach. Eor 
this reason we need to perform hydro dynamical simula¬ 
tions in order to provide a better description of the main 
physical events in a nuclear region where a SB coexists 
with an AGN. 

4. THE NUMERICAL MODEL 
4.1. Seyfert galaxy initial setup 

In our model we consider the central (kpc) region of 
a typical spiral (Seyfert) galaxy and set the initial con- 
ditions for the ISM following t he procedure outlined in 
IMelioli et al.l (j2QQ8l . I2QQ9L l2Q13f ) . We first assume a mass 
model for the galaxy, which includes the contribution of 
a bulge and a stellar disk, and then we set the ISM in 
equilibrium with the gravitational potential given by the 
summation of the dark matter halo, the bulge, and the 
disk contributions. The dark matter halo is assumed to 
follow the Navarro, Erenk & White profile (iNavarro et al.l 
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Log SFR,^,(Mg yr-‘) 


Figure 3. Mass loss rate for a nuclear region of a galaxy consid¬ 
ering an active volume with a radius R=300 pc and a thickness h= 
200 pc as a function of the SFR, for different values of the heating 
efficiency, HE. Solid (black) line: HE = 3%; long-dashed (red) line: 
HE = 10%; dashed (blue) line: HE = 50%; dot-dashed (magenta) 
line: HE = 100%. The SFR and the mass loss rate are given in 
M© yr“^ (logio scale). 


^dm(r) = - 


GMvir ln(l + x) 

Fs/(c) X 


( 8 ) 


where r is the spherical radius, Mvir is the mass at the 
virial radius, Tvir defined as the radius where the aver¬ 
age density is ^ 10^ times the cosmological critical den¬ 
sity, Pcrit, Fs is a scale radius, x is the ratio between 
r and rg, c = rvir/^s is the concentration and /(c) = 
ln(l + c) — c/(l -h c). In this model we considered Mvir 
= 10^^ M©, Tvir = 120 kpc, Tg = 10 kpc and c = 12, 
adopting a A cold dark matter cosmological universe with 
(1 m= 0-27, r2A=0.73 and Ho=71 km s“^ Mpc“^. How¬ 
ever, we note that due to the limited extension of the 
galaxy region here investigated, we verify that the gravi¬ 
tational potential of the dark matter halo has almost no 
influence in the determination of the initial distribution 
and the evolution of the gas. 

The gravitational potenti al of the stellar b ulge is given 
by a Plummer distribution (jPlummerlllQllf) : 


^h{r) 


GMi, 

Vc + r 


(9) 


where Mb is the total mass of the bulge and Vc is the 
core radius defined as the radius containing 50% of the 
stellar mass. The bulge mass is assumed to be always 
Mb = 10^^ M©, while the core radius varies according to 
the model. 

Finally, the gravitational potential of the gas 
disk is given by a M iyamoto & Nagai profile 
(|Mivamoto fc Na^ll975f ). that is: 


^disk{r,z) 


GMdisk 

r‘^ [az‘^ -fb^) ^ 


( 10 ) 


where Mdisk is the mass of the disk (for the whole galaxy), 
while a and b are its radial and vertical scales, respec¬ 
tively. In this study we assume Mdisk = 5 x 10^^ M©, a = 
8 kpc and 6 = 0.8 kpc. Therefore the total gravitational 
potential is: ^{r,z) = <I>dm + ^b{r,z) -b ^disk{r,z). 

The rotating gas in the disk is initially put in equi¬ 
librium with the galaxy gravitational potential. We 
obtain the rotation velocity resolving the equation 
v(r,z) = y(r di » (r,2:)) /dr)z=r.onstant, fsee iMelioli et alj 
12008112009L1201311201511 . The values of the velocity ob- 
tained in this way mimic quite well the rotation curve of 
a spiral galaxy. 

We assume three distinct galaxy models, SyL, SyM 
and SyH characterized by gas disk distributions with 
different total mass, density profile and resulting column 
densities defined as low {Nh = 10^^ cm“^), medium 
{Nh = 10^^ cm“^) and high {Nh = 10^^ cm“^), respec¬ 
tively. All models have a multi-phase gas distribution, 
and each phase at t = 0 is defined by a typical tempera¬ 
ture Ti^Q and central density initially in equilibrium 
in the total gravitational potential described above. The 
isothermal density distribution of each phase will have 
the form: 


Pi{r,z) = po^i exp 


$(r, z) 


e|$(r, 0) - (1 - e?)$(0,0) 

disk,i 


( 11 ) 

where ■ is the isothermal sound speed of the Tphase 
of the gas and quantifies the fraction of rotational sup¬ 
port of the ISM. In our model the i-phase density is re¬ 
placed by the (i + l)-phase density wherever the {i + 1)- 
phase pressure is larger than the (coldest and densest) 
pressure of the i-phase. 

The galaxy parameters adopted in this study are pre¬ 
sented in Table [H and the initial density and temperature 
distribution for each model is shown in Fig. HI 


4.2. Energy injection 

As stressed before, the energy injected in the central 
region of a Seyfert galaxy comes from a large number of 
sources: stellar winds and protostellar jets, SN (type I 
and II) explosions, SMBHs, cosmic rays and spiral waves. 
However, the average energy injected by the protostellar 
jets is about 100 time s smaller than the energy i njected 
by the SN explosions (|Mac Low fc Klessenll2QQ4 ). while 
only the Wolf-Ra yet stars have winds th at are energet¬ 
ically important (|Nugis fc LamersI |2000[) . but last only 
^ 10^ years. Furthermore, since we are mainly inter¬ 
ested in the gas outflow, cosmic rays and spiral waves 
may be also neglected because although they may play 
an important role in the star formation processes, they 
are expected to be less relevant to drive the gas out of 
the disk in active galaxiefl 

^ See, however, the importance of cosmic rays to drive winds in 
normal spiral galaxies and dwarf galaxies (see, e.g.. lEverett et ^ 
[200^ : [B5oth et ahllMT^l 
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Table 1 

Parameters for the galaxy models setup SyL (thin disk), SyM (intermediate disk) and SyH (thick disk) 


Model 

Tc 

kpc 

P0,1 

cm“^ 

P0,2 

cm“^ 

P0,3 

cm“^ 

Lzisfe, 1 
K 

LizsA:,2 

K 

LizsA:,3 

K 

ei 

62 

63 

Nh 

cm“^ 

Maoopc 
10® M0 

SyL 

1.3 

1 

10-^ 

10-=^ 

10-^ 

5 X 10^ 

10® 

1 

0.9 

0.8 

2x10^1 

0.03 

SyM 

0.7 

30 

1 

3x10-2 

10“^ 

5 X 10^ 

10® 

1 

0.9 

0.5 

4x10^2 

0.5 

SyH 

1.3 

100 

1 


10^ 

5 X 10^ 

10® 

1 

0.9 

0.5 

2 X1023 

3.0 


SyL 

-0.00 
- 0.66 
- 1.31 
- 1.97 
- 2.62 
- 3.28 
- 3.94 

- 0.4 - 0.2 0.0 0.2 0.4 

R 

SyM 


0.4 

0.3 

M 

0.2 

0.1 


- 0.4 - 0.2 0.0 0.2 0.4 

R 





Figure 4. Edge-on logarithmic density distribution (with density 
given in cm“^) of the unperturbed ISM for the models SyL (top 
panel), SyM (middle panel) and SyH (bottom panel). Distances 
are given in kpc. 

For these reasons, in this study we assume that the 
kpc-scale outflow observed in different wavelengths is 
mainly driven by the Type I and Type II SN explosions 
and/or by the jet produced by the activity of the SMBH 
in the center of the galaxy. 

The rate of type I SN explosions is proportional to the 
bulge mass (M^), which in turn is proportional to the 
mass of the central SMBH (Mbh) (see §3.1). Given th e 
relation Mbh/Mi) ^ 10“^ (see, e.g.. lHaring L Rixll2QQ4l ). 
and assu ming a SNI rate in the bulge of ^ 0.01 yr“^ per 
10^^ Mq (|Pain et all 1 1 9961) . we obtain: 


random time ti in the interval 0 < < tsF, where tsF=30 

Myr is the time of stellar activity. This procedure gener¬ 
ates points with a random spatial frequency proportional 
to the stellar de nsity described by the Plummer profile 
Inmmer 


Pb,star('^) — 


3Mb 

47rrc 


-5/2 


(13) 


To introduce the SNII (characterized by a lifetime < 
30 Myr), we assume a SB in the central region of the 
host galaxy, in the disk within a radius of ^ 300 pc and 
with a typical half-height (above and below the disk) 
of 80 or 160 pc, depending on the model. We con¬ 
sider an event of star formation in which a total mass 
of stars ^ 3 X 10^ or 3 x lO^M© (depending on the 
model) form. Considering a SN rate is ^ 10“^ (or 10“^) 
yr“^, which is in agreem ent with many observa tions of 
Seyfert galaxies (see, e.g.. lForbes fc Norrislll998[) . we ex¬ 
pect to have ^ 3 xlO^ (or 3 xlO^) SNe injecting an 
average luminosity Lsn ^ 3x10^^ (or 3x10^^) erg s“^ in 
a time interval of 30 Myr. As in the case of the SNI, to 
set the SNII spatially and temporally within the SB re¬ 
gion, we associate randomly to each /-th SNII a position 
P* = (ri,hi), where 0 < < 300 pc and (-80,-160) 

pc < hi < (80,160) pc depending on the model, and a 
time in the range 0 < < 30 Myr. Both for the SNe 

I and the SNe II, we inject the energy of each explosion 
over a tir ne of 300 yr and with a SN heating effi ciency HE 
= 100% (jMelioli fc de Gouveia Dal Pinoll2003 ) FI. There¬ 
fore, each SN explosion injects mass and energy at rates 
M = Minj/(300 yr) and Lsn = 10^^/ (300 yr), respec¬ 
tively, where Minj is the mean mass released by a single 
SN explosion, that is, 16 Mq for the SNII and 1 Mq for 
the SNI. 

Finally, we assume that the SMBH injects a total con¬ 
stant luminosity of ^ 10^^ erg s“^ above and below the 
midplane of the galaxy in two collimated jets character¬ 
ized in most of the models by a non-relativistic velocity 
perpendicular to the plane of the disk of 2.1 x 10^ cm s“^ 
and by a total rate of injected matter of ^ 6 xlO ^ Mq 
yr“^. The jet is injected in a central cylindric volume of 
about 14 pc^. In order to verify how the jet setup may 
affect the results, we have also considered a model with 
a higher velocity and smaller density jet, with the same 
luminosity of the one described above, but with a rela¬ 
tivistic velocity of ^ 0.2 c and a rate of injected matter 
of ^ 5 xl0“^ Mq yr“^. 

Since in this study we are interested in understanding 
the role of each of the energy sources described above on 


SNI = 10-^Mbh yr“^ (12) 

where the mass of the SMBH, Mbh^ is given in units of 
Mq. Each i-th SN explosion is associated randomly to a 
given position Pi in the stellar bulge of the galaxy at a 


We note that a SNR will develop only after the SN shock front 
enter the Sedov phase, that is, after a time ~ 400/n^/^ yr. The 
SNR evolution may be considered adiabatic up to a time tgedov ~ 
3 X 10^yr. This explains why we adopt HE = 100% in the 
first 300 yr. 




















the galaxy evolution we have performed different simula¬ 
tions where either all or only part of the energy sources 
have been considered. In order to cover a parametric 
space as large as possible, we have run 16 different mod¬ 
els and a summary of the characteristics and parameters 
adopted for each model are given in Table O 

4.3. Numerical methodology 


To simulate the evolution of the system descrived 
above, we employed a modified version of the Carte¬ 
sian Godun ov MHD code originally developed by G. 


Kowal fsee iKowal et al. 

120071; 

lEalceta-Goncalves et al.l 

I 2 OO 8 I: iSantos-Lima et al. 

I 2 OIOI: iLeao et al.l 2013D. nsiner 


the Harten-Lax-van-Leer C (HLLC) solver and a sec¬ 
ond order RungeKutta for time integration. It also uses 
the message-passing interface (MPI) library to achieve 
portability and efficient scalability on a variety of paral¬ 
lel high-performance computing systems. The code is 
available upon request directly from the authors. 

In this version, the code includes a parametrized cool¬ 
ing function in the energy equation that allows the gas 
to cool down to 10 K with errors smaller than 10% and 
which is calculated implicitly in each time step for each 
grid position. The cooling function considers an opti¬ 
cally thin gas in ionization equilibrium and takes into 
account also the gas metallicity and the gas fraction of 
H 2 , a c cordi ng to the methodology used by iRaga et~all 
(I 2000112001 . We have run all the models with a maxi¬ 
mum resolution of 1.9 pc per cell. 0 The adopted box 
has physical dimensions of 1x1x1 kpc in the x, y and z 
directions, respectively, and is covered by 512^ cells. 

The energy associated to each SN explosion is injected 
as thermal energy in a single cell, while the SMBH jet 
is injected in a single cell above and below the midplane 
disk, at the center of the system. 

In all simulations the initial metallicity at the mid¬ 
plane is assumed to be equal to the solar one {z/ zq=1) 
and decreases inversely proportional with the height of 
the disk, up to a minimum value Zmin /We note 
that the metallicity can influence the value of the cooling 
function up to a factor 10 specially at gas temperatures 
^ 10^ K, but it is not so important for the gas cooling 
at temperatures between 10^ and 5 x 10^, or above 10^ 
K. For this reason, we have adopted for the injected SNe 
a metallicity that inc reases accord i ng to their evolution 
as described, e.g., in iMelioli et al.l (|2Q13L see also refer¬ 
ences therein). We further notice that we have also run 
a few models (not presented here) with different initial 
values of the metallicity in the disk (z/2;©=0.3; zjz^—N) 
and found no significant modifications in the results with 
regard to the reference model {z/Z q=1). 

Each simulation was run on 512 processors of the LAI 
supercomputing centre at lAG-USP, using about 10^ 
GPU hours per simulation. 

5. RESULTS 

In this section, we present the main results for the 
different models listed in Table 2. All the models were 

^ We have also performed several tests with a larger resolution of 
0.95 pc per cell for about 2 Myr and found that the results are very 
similar to those with a 1.9 pc resolution per cell. For this reason, 
we adopted the latter scheme for the study presented here as it 
allowed us to save computation time without loosing information 
and quality in the results. 


followed for a time interval between 2 and 12 Myr (de¬ 
pending on the model), which is an appropriate time to 
investigate the global properties of the wind build-up, 
from the first SN explosions up to a nearly steady-state 
phase, when no gas inflow (from the external regions to 
the center of the galaxy) is detected. 

Although we analysed all models of Table 2, we will 
focus our attention on four of them, namely, the models 
SyM-SNI-SB, SyM-SNI-SB-JET, SyH-SNI-SB and SyH- 
SNI-SB-JET which combine one or more of the poten- 
tial outflow-dr i ving m echanisms. Observations (see, e.g., 
iTombesi et al.l 12013D indicate that a large fraction of 
Seyfert galaxies has an amount of gas in the inner region 
of the order of 10^ M©, so that models with the SyH 
setup (Table 1) can be taken as main references in this 
study. On the other hand, models with the SyL setup, 
which have lower values of gas column density in the 
center, are more rare to observe and for this reason their 
results will be only briefly discussed in the Appendix. 

5.1. Seyfert galaxy models with intermediate column 
density (SyM) 

5.1.1. SyM-SNI 

In this model, we consider a disk galaxy with the SyM 
setup (see Table 1), characterized by a total mass in the 
central half kpc^— volume of about 5 xlO^ M©, a max¬ 
imum density of 30 cm“^ and a column density, nn = 
10^^ cm“^. The energy injected comes only from the SNI 
explosions, which are expected to occur in any galaxy re¬ 
gardless of the presence or not of a SB region or a SMBH 
jet. In this case we find that the SNe I, alone, are unable 
to remove the gas from the center of the galaxy. Eigure[6] 
depicts the mass evolution for this model (dashed-green 
line) and we clearly see in the top panel that the mass 
of the whole system remains nearly constant over the 
time simulated, that is, 6.5 Myr for this mode0 In the 
core of the galaxy, within a radius of 40 pc, the mass of 
the gas varies by about ± 15% compared to its initial 
value, while in the disk up to a radius of 400 pc the total 
amount of gas is reduced by only a few per cent (see Eig. 
[6l dashed-green lines in the middle and bottom panels). 

Eigures [5] and [8] (dashed-green line) indicate that the 
SNI explosions only induce the growth of turbulence in 
the disk and the halo. However, there is no impact in 
the evolution of the thermal pressure or the maximum 
density (which are almost constant along the simulation), 
so that no star formation process, gas outflow or inflow 
occurs in this case (see also Eig. [7]). We note that in 
Eigure [51 the average velocity refers to the whole box 
and this average (performed over the volume and not 
over the mass) is strongly affected by the high velocities 
in the halo, above and below the disk. However, the halo 
(due to its very low density) does not contribute to the 
evolution of the disk gas and for this reason, despite the 
high average velocity, the system does not generate any 
significant outflow. 

® We note that in this and all diagrams representing the global 
evolution of the system, there is a small decrease in the gas mass 
mainly due to the rotation of the galaxy which causes an artificial 
mass loss through the lateral boundaries of the computational do¬ 
main. For this reason, in our simulations, in order to compute the 
final mass loss due to the stellar (or jet) feedback we calculate the 
difference between the total mass in a model with gas outflow and 
its counterpart without gas outflow. 
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Table 2 

Initial conditions adopted in the hydrodynamic simulations for 16 different models. For each Model: a) 
Name; b) Host galaxy setup; c) SNI rate (yr“^); d) SF rate (Mq yr“^); e) SB region half-hight (pc); /) 
SMBH Jet luminosity (erg s“^); g) SMBH Jet velocity (vjet/c) 


Model® 

Gal. Type** 

SNF 

SFR<i 

ue 

^SB 

Jet Lum.t 

Jet Vel.^ 

SyL-SNI 

SyL 

0.01 

- 


- 

- 

SyL-SNI-SB 

SyL 

0.01 

1 

160 

- 

- 

SyL-SNI-SB-JET 

SyL 

0.01 

1 

160 

10^2 

5x10-2 

SyL-SNI-JET 

SyL 

0.01 

- 

- 

1042 

5x10-2 

SyM-SNI 

SyM 

0.01 

- 

- 

- 

- 

SyM-SNI-SB 

SyM 

0.01 

1 

80 

- 

- 

SyM-SNI-SB-JET 

SyM 

0.01 

1 

80 

(M 

o 

1—1 

5x10-2 

SyM-SNI-JET 

SyM 

0.01 

- 

- 

1042 

5x10-2 

SyM-SNI-lOSB 

SyM 

0.01 

10 

80 

- 

- 

SyH-SNI 

SyH 

0.01 

- 

- 

- 

- 

SyH-SNI-SB 

SyH 

0.01 

1 

80 

- 

- 

SyH-SNI-SB-large 

SyH 

0.01 

1 

160 

- 

- 

SyH-SNI-SB-JET 

SyH 

0.01 

1 

80 

(M 

o 

1—1 

5x10-2 

SyH-SNI-JET 

SyH 

0.01 

- 

- 

1042 

5x10-2 

SyH-SNI-SB-JET-light 

SyH 

0.01 

1 

160 

- 

2x10-1 

SyH-SNI-lOSB 

SyH 

0.01 

10 

80 

- 

- 


5.1.2. SyM-SNI-SB 

Also in this model the disk galaxy is built with the 
SyM setup (see Table 1), but the energy sources now 
include both regular SNI explosions and those of a SB 
region characterized by a SFR of 1 Mq yr“^. We calcu¬ 
late the disk evolution over a time of 12 Myr (Fig. [9]) and 
find that between 2 and 6 Myr an average gas outflow of 
about 15 Mq yr“^ is driven when the SB region is active 
(see Figures [6] and 0 dashed blue lines). The mass of the 
disk decreases by about 4.5 x 10^ Mq in the first 5 Myr 
and then, since there is no gas infall coming from the ex¬ 
ternal region of the galaxy (not considered in this study) 
the gas density drops to values of about 10“^ cm“^ (see 
Fig. [H dashed blue lines), and a fraction of the gas raised 
by the SN explosions falls back to the disk, leading to a 
final average density, after 12 Myr, of about 1 cm“^ (Fig. 
[ 91 ). After 2 Myr the gas in the core of the galaxy (within 
a radius of 40 pc) is completely removed, and globally 
the kpc^ volume here considered looses about 3.5 x 10^ 
Mq of gas along 2 Myr, corresponding to a total mass 
transfer rate of 1.75 Mq yr“^. 

Observing Fig. [9] we also see that between the hot 
low density gas outflow and the highly energized gas in 
the disk, there is a third warm phase characterized by 
a temperature of 10^ K, the minimum allowed in our 
model. This gas, which enters in the hot wind and flows 
out of the system, far away from the ionizing sources of 
the galaxy may be able to cool to low temperatures and 
eventually, in presence of dust, may form molecular gas 
at few kpc above t he disk, as it has be en observed in 
some systems (e.g., iMorganti et al.ll2Q13L ; see also the 
discussion section 6.2). 

5.1.3. SyM-SNI-SB-JET 

This model is equal to the previous one (SyM-SNI- 
SB), but in this case we considered also the presence 
of a highly collimated jet (Table 1) injected in a cylin¬ 
drical volume characterized by a height of 1.9 pc above 
and below the midplane of the disk, a mass at a rate 
of about 6 X 10“^ Mq yr“^, a (non-relativistic) velocity 
of about 2x10^ cm s“^ and a total kinetic luminosity 
of 10^^ erg s“^. We calculated the evolution of the sys¬ 


tem over a time of 6 Myr. The results, shown in Figs. 
ElE and [8] (solid-black lines) and Fig. [TOl demonstrate 
that the overall gas evolution is basically the same ob¬ 
tained without the jet. The only significant difference 
is in the central core of the disk, within 40 pc, where 
the gas mass is removed faster, in about 10^ yr. The 
jet at the beginning of its evolution generates a trans¬ 
verse shock wave that propagates radially until reaching 
an equilibrium position at about 30 pc from the galaxy 
center. This makes the gas in this region denser and less 
vulnerable to the action of the SN explosions, which re¬ 
duces and delays the total gas outflow. The presence of 
the jet, as expected, also causes a higher mean velocity 
in the whole system and a very high velocity signature in 
the low density gas above the disk. However, the mass, 
density, pressure and temperature evolution of the disk 
are comparable with those of the model SyM-SNI-SB, 
and therefore we may conclude that in this case the jet 
does not affect significantly the global evolution of the 
system. 

5.1.4. SyM-SNI-JET 

In this model we consider a disk with a SyM setup, 
the presence of SNI explosions in the bulge and a central 
highly collimated jet with a total kinetic luminosity of 
10"^^ erg s“^. Differently from model SyM-SNI-SB-JET, 
in this case we do not consider a SB region. The results of 
the evolution of this model over 6 Myr (Figs. |6l0and|Hl 
short-dashed purple lines) show, as in the previous case, 
that the jet alone is unable to remove from the system 
a significant amount of gas. Indeed, the fraction of gas 
shocked by the jet along the axis of propagation corre¬ 
sponds to about 2% of the gas in the nuclear region of the 
galaxy (within R < 300 pc), but there is a higher increase 
in density and pressure in the galactic disk induced by 
the jet in its normal direction. As in the previous model, 
during the first 2 Myr only, the gas in the core of the 
galaxy (within R < 40 pc) is swept away by the shock 
wave generated at the beginning of the jet propagation. 
However, this phenomenon is local and transient, and 
does not affect much the global evolution of the system. 
After 2 Myr, about 30% of the gas falls back into the 
nuclear region, and the evolution of the galaxy continues 
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Figure 5. Edge-on logarithmic gas density (top panel) and ve¬ 
locity (absolute value, bottom panel) distribution at t = 4.2 Myr, 
for the model SyM — SNI. Distances are given in kpc, density 
is in cm“^ and velocity is in units of the reference sound speed 
computed at T=5xl0^ K, 5 x 10 ^ = 33 km s“^. 

to be the same as in Model SyM-SNL 

5.2. Seyfert galaxy models with high column density 

(SyH) 

5.2.1. SyH-SNI 

In this set of models we consider a disk galaxy with 
the SyH setup (see Table 1), characterized by a total 
mass in the central half kpc^ of about 3 xlO^ Mq, a 
maximum density of 100 cm“^ and a column density nn 
= 2 xlO^^ cm“^. In this configuration the disk is thicker 
and more massive than in the previous setup, resembling 
the gas distribu tion observed in th e nuclear region of 
Seyfert galaxies (jTombesi et alJl2013[) and also of our own 
Galaxy (see, e.g., Vergani et ah 2004). 

As in model SyM-SNI, in the model SyH-SNI the only 
energy source comes from the SNI explosions that alone 
are unable to remove the gas from the center of the galaxy 
(Fig. E}. The mass of the system is nearly constant over 
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Figure 6. Time evolution of the mass of the gas within the whole 
system (z < 500 pc, upper panel), the thick disk (z < 200 pc, 
middle panel) and the central core of the galaxy (r < 40 pc, bottom 
panel) for the models SyM-SNI (long dashed-red lines), SyM-SNI- 
SB (dashed-blue lines), SyM-SNI-JET (dot-dashed-magenta lines) 
and SyM-SNI-SB-JET (solid-black lines). Time is in Myr and mass 
is in units of Mq, logarithmic scale. 
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Figure 7. Time evolution of the gas mass transfer and loss rate of 
the thick disk (z< 200 pc) for the models SyM-SNI (long dashed- 
magenta lines), SyM-SNI-SB (dashed-blue lines), SyM-SNI-JET 
(dot-dashed-red lines) and SyM-SNI-SB-JET (solid-black lines). 
Time is in Myr and mass loss rate is in units of Mq yr“^. 
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Figure 8. Time evolution of the main physical variables within 
the system: maximum number density (upper panel), average ve¬ 
locity (middle panel) and average pressure (bottom panel). As 
in Figs. [ 6 ] and [Jl we have the models SyM-SNI (long dashed- 
red lines), SyM-SNI-SB (dashed-blue lines), SyM-SNI-JET (dot- 
dashed-magenta lines), and SyM-SNI-SB-JET (solid-black lines). 
Density and pressure are in cgs, velocities are in km s“^ (all vari¬ 
ables are in log-scale). 

the maximum simulated time (9 Myr) aud uo siguificaut 
mass trausfer from the core of the galaxy to above the 
disk {h > 200 pc) is detected, as depicted iu Figures [12] 
aud Fig. [131 Therefore, as expected, the thick disk iu a 
uuclear regiou of a galaxy, eveu if completely photoiou- 
ized by the ceutral source (a SMBH iu the case) aud 
thus characterized by a miuimum temperature of 10^ K, 
is uot destroyed or eveu partially removed by the SNI ex- 
plosious, which however coutribute for maiutaiuiug the 
turbuleuce of the ISM, maiuly above the disk (Fig. [T4|) . 

5.2.2. SyH-SNI-SB 

Iu this case the euergy iujected iu the thick disk {SyH) 
comes from the SNI aud SNII explosious associated to 
the bulge aud to a SB regiou with a SFR of 1 Mq yr“^, 
respectively. Siuce the thick disk exteuds up to a height 
hz ^ 180 pc, the stars withiu the SB regiou may be dis¬ 
tributed iu differeut ways. Iu this study we cousider two 
differeut stellar distributious iu a uuclear regiou with a 
radius Rsb = 300 pc. Iu oue case they are uuiformly 
distributed betweeu a height -80 < < 80 pc, while 

iu the secoud case (here defiued large) the stars are dis¬ 
tributed with a power law oc {\hz\/ho)^ betweeu a height 
-160 > hsB ^ 160 pc aud with ho = 160 pc, that is, 
the SB regiou is active throughout the thickuess of the 
disk. We have calculated the disk evolutiou over a time 
of 9 Myr aud we see that the global evolutiou is quite 
similar iu both cases (Figs. [T5] aud [T6|) . This meaus that 
the SNII vertical distributiou, withiu the disk, is uot so 
importaut as the total SNII uumber aud/or the SNII 
deusity surface, aud for this reasou iu the uext models 
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Figure 9. Edge-on logarithmic gas density (npper panel) and 
velocity (absolnte valne, vertical component, bottom panel) distri- 
bntion for the model SyM-SNI-SB at t = 5 Myr. Density in cm“^ 
and velocity in nnits of the reference sonnd speed (c^ 5 x 10 ^ = 33 
km s“^). 


we will cousider ouly the large stellar distributiou withiu 
the SB which is possibly a more realistic case. 

Iu both models above, after about 2.5—3 Myr, the gas 
euergized by the SN explosious begius to rise from the 
disk to the halo at au average rate of ^ 20 Mq yr“^, aud a 
maximum rate of 45 Mq yr“^ at t ^ 4 Myr. Aloug 7 Myr 
the mass of the disk iu the large case decreases by about 
10^ Mq, that is, ^ 30% of its iuitial mass is lifted above 
the deuser disk to heights > 200 pc. Similarly, wheu the 
SNII are uuiformly distributed withiu a thiuuer regiou, 
the mass of the disk is reduced by ^ 23%, with au active 
gas outflow startiug at t = 3 Myr aud eudiug after 3 
Myr (i.e., at t = 6 Myr). Also the deusity, the pressure 
aud the velocity evolutious are very similar iu both cases 
(Fig. [TH The ouly differeuce is observed iu the removal 
of gas from the uuclear regiou, clearly more efhcieut iu 
model SyH-SNI-SB. Iu this case, all the gas is dispersed 
to larger radius, while iu the large case about 10% of 
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Figure 10. Edge-on logarithmic gas density (top panel) and ve¬ 
locity (absolute value, vertical component, bottom panel) distribu¬ 
tion for the model SyM-SNI-SB-JET at t = 5 Myr. Density is in 
cm“^ and velocity is in units of the reference sound speed (c^ sxio^ 
= 33 km s“^) 

the gas remains in the central 40 pc. However, these 
differences are not enough to affect the evolution of the 
system and we can state that the choice of the vertical 
distribution of the SNs does not seem to be so important 
in this study. Close to the midplane of the disk, the 
SB activity generates a 2-phase ISM characterized by a 
diffuse gas with density of ^ 10“^ — 10“^ cm“^ and 
temperature of about 10^ K, where dense clumps and 
filaments are embedded with a density of ^ 10 cm“^ 
and temperature of 10^ K. 

Globally, the mass of the system here considered (with 
a kpc^ volume) decreases by about 10^ Mq at an average 
rate of few Mq per year (top panel of Fig. [12]), but since 
the region we have simulated extends only up to ± 500 
pc (on each side of the disk), we are unable to predict 
the final fate of this gas, i.e., whether it will fall back 
to the disk or definitively leave the galaxy as a wind. 
An estimate of the inner galactic disk escape velocity 


Figure 11. Edge-on logarithmic gas density (top panel) and gas 
velocity (bottom panel) distribution at t = 4.2 Myr, for the model 
SyH-SNI. Distances are given in kpc, density is in cm“^ and veloc¬ 
ity is in units of the reference sound speed computed at T=5xl0^ 
K, Cg 5 xio4 = 33 km s“^. We note that here the velocity signal 
characterizes the real direction of the velocity (rather than its ab¬ 
solute value). It is generally positive above the disk and negative 
below it, characterizing an outflow motion. We also identify a few 
patches where the velocity has an opposite signal with regard to 
the surrounding gas, characterizing regions of fallback flow. 

((2GM/r)^/^) indicates that it is of the same order of 
the outflow velocity of part of the energized gas, so that 
both fates are possible. 

5.2.3. SyH-SNI-SB-JET 

As in the study of the SyM setup, this model has sim¬ 
ilar initial conditions to the previous one (SyH-SNI-SB), 
but includes also the presence of a highly collimated jet 
with the same characteristics described above (Table 1 
and §4.2 4.2). We calculated the evolution of the system 
over a time of 6 Myr and the results are shown in Figures 
[121 [131 (HI (solid-black lines), and [T71 Also in this case, 
the gas evolution does not differ much from the model 
above without the jet. The gas transport above the disk 
is somewhat less efficient than in the SyH-SNI-SB large 
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Figure 12. Time evolution of the mass of the gas in the whole 
system (z < 500 pc, upper panel), in the thick disk (z < 200 pc, 
middle panel) and in the central core of the galaxy (r < 40 pc, bot¬ 
tom panel) for the models SyH-SNI (long dashed-red lines), SyH- 
SNI-JET (short-dashed-magenta lines), SyH-SNI-SB-JET (solid- 
black lines), SyH-SNI-SB large (dashed-blue lines), and SyH-SNI- 
SB (dot-dashed-purple lines). Time is in Myr and the mass is in 
units of Mq, logarithmic scale. We note that the earlier truncation 
in some tests was due to computational time constraints. From ear¬ 
lier tests with lower resolution we have verified that the results did 
not change substantially, then we truncated the longer lasting runs 
earlier. 

model, with a mean mass transfer rate of about 15 Mq 
yr“^ (and a maximum of ^ 30 Mq yr“^), while in the 
central core of the disk, within 40 pc, about 90% of the 
gas is removed during the first 4x 10^ yr. due to the pres¬ 
ence of the jet. However, this trend does not represent a 
steady state and we note that after 2 Myr the gas evolu¬ 
tion is very similar to that of the previous model. Also 
the density, pressure and temperature evolution of the 
disk are comparable with those of the model SyH-SNI- 
SB large, so that we may conclude that also in a system 
with properties that resemble those of a Seyfert galaxy 
with a thick disk, the SMBH jet is not expected to affect 
the global evolution of the gas in the inner regions, up 
to r\j 500 pc. 

To understand better the evolution of this system SyH- 
SNI-SB-JET compared to the previous one without the 
jet (SyH-SNI-SB large), we have plotted in Figs. [Island 
[19] the time evolution of the gas mass for different values 
of the vertical velocity and the temperature, respectively. 
We note that the gas mass related to a given velocity 
and temperature in the model including the jet repre¬ 
sents only a few per cent of the total mass of the nuclear 
region. This explains why we see no significant differ¬ 
ences in the gas evolution between the two models. In 
other words, though the presence of the jet may change 
a little the velocity and temperature profiles and thus 
the emission structure, it does not seem to be efficient 
enough to remove much gas at the galactic scales. 
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Figure 13. Time evolution of the gas mass transfer and loss rate 
of the thick disk (z< 200 pc) for the models SyH-SNI (long-dashed- 
red lines), SyH-SNI-JET (short-dashed-magenta lines), SyH-SNI- 
SB-JET (solid-black lines), SyH-SNI-SB large (dashed-blue lines), 
and SyH-SNI-SB (dot-dashed-purple lines). Time is in Myr and 
mass loss rate is in units of Mq yr“^. 
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Figure 14. Time evolution of the main physical variables within 
the system: maximum number density (upper panel), average 
velocity (middle panel) and average pressure (bottom panel) for 
the models SyH-SNI (long-dashed-red lines), SyH-SNI-JET (short- 
dashed-magenta lines), SyH-SNI-SB-JET (solid-black lines), SyH- 
SNI-SB large (dashed-blue lines), and SyH-SNI-SB (dot-dashed- 
purple lines). Density and pressure are in cgs, velocities are in km 
s“^ (all variables are in log-scale). 
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Figure 15. Edge-on logarithmic gas density (upper left panel), temperature (upper right panel), pressure (bottom left panel ) and velocity 
(vertical component, bottom right panel) distribution for the model SyH-SNI-SB large at t = 3.75 Myr. Density is in cm“^, temperature 
in K, pressure in units of 7x10“^^ erg cm“^ (code units) and velocity in units of the reference sound speed (c^ sxio^ = 33 km s“^). Here 
the velocity signal characterizes the real direction of the flow velocity. 


5.2.4. SyH-SNI-JET 

To further stress the results and conclusions above and 
verify whether or not in the absence of a SB the remain¬ 
ing injecting mechanisms are still able to sweep large 
amounts of gas from the thick disk to the halo, we have 
also run a model including only the jet and the SNL The 
magenta short-dashed lines of Figs. [121 [l3] and [14] indi¬ 
cate that the only effect of the jet on the evolution of 
the system occurs in the nuclear region, within 40 pc, 
as in the previous model. No gas mass is transported to 
outside and no gas outflow from the disk to the halo is 
observed. The jet crosses the disk affecting it just in the 
very beginning of its evolution. After a few hundreds of 
thousands years the jet interaction with the disk ceases 
and it continues to flow freely through the low-density 
channel previously opened by it. 

5.2.5. SyH-SNI-SB-JET-light 


As remarked in §4.2, aiming at verifying how the char¬ 
acteristics of the jet can affect the evolution, we have also 
considered a lighter jet with a relativistic velocity of ^ 
0.2 c and a rate of injected matter of ^ 5 xl0“^ Mq 
yr“^. In this case, though the jet speed, Uj, is larger, 
the density contrast between the jet and the disk near 
the base is much smaller (77 = pyet/Pdis/c=5.5x 10 “^) 
than in the previous model (SYH-SNI-SB-JET). Consid¬ 
ering the momentum flux conservation at the shock that 
forms at the head of the jet, one can easily demonstrate 
that the jet propagation velocity into the disk, de¬ 
pends both on the jet speed and the density contrast, i.e. 
Vh — (e.g. Ide Gouveia dal Pino fc Benz! 

Il993[ ) . This indicates that the larger the density contrast 
and the jet speed the faster the jet propagates. For the 
model above with a light jet (SyH-SNI-SB-JET light) ^ 
the much smaller density contrast gives a propagation 
velocity vy smaller than that of model SyH-SNI-SB-JET, 
in agreement with the results of the simulations (see Fig. 
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loss rate of the models with (SyH-SNI-SB-JET light) and 
without (SyH-SNI-JET light) a SB region. 

Eurthermore, since the faster the jet crosses the disk 
less energy is deposited into the ISM of the host galaxy, 
the jet with the higher density contrast (SyH-SNI-SB- 
JET in Eigure [20l right panels) produces a more colli¬ 
mated flow. On the other hand, the longer time of energy 
deposition into the disk by the lighter jet in SyH-SNI- 
SB-JET light model due to the slower evolution, results 
a similar mass evolution to that of Model SyH-SNI-SB- 
JET, with differences of only a few Mq per year for the 
mass transfer rate of the thick disk and no significant 
differences in the mass loss rate of the whole system. 

5.3. Seyfert galaxy models with higher SFR (lOSB) 

Einally, we have run models (not shown here) con¬ 
sidering a SER of the SB region 10 times higher than 
in the previous cases. We considered a SyH setup 
with (SyH-SNI-lOSB-JET model) and without (SyH- 
SNI-IOSB model) the jet. These models were run either 
with higher or lower gas metallicity than the solar val- 
uefl and with different SN distribution, considering both 
single and clustered SN explosions. We have found simi¬ 
lar results between these tests, nearly independent of the 
metallicity (see, hoever, iMelioli et al.ll2QQ9l I2Q13LI2Q15L 
for a more detail study on the effects of metalicity in 
galactic and SB winds), which are also consistent with 
those presented above. Of course, the mass transfer rate 
from the disk to the halo is higher and, almost insensi¬ 
tive to the presence of the jet, reaches a maximum value 
of ^ 170 Mq yr“^ and a mean value of ^ 90 Mq yr“^, 
resulting in a total mass exchange between the disk and 
the halo involving about 40% of the mass of the system. 
Moreover, the higher energy injected by the SB region 
increases the thermal pressure of the gas in the nuclear 
region and consequently, the channel dug by the jet is 
more confined and tighter. This result suggests that the 
higher the SER and the column density of the disk, the 
greater is the collimation of the SMBH jet. 


Figure 16. Edge-on logarithmic gas density (top panel), and 
velocity (vertical component, bottom panel) distribution for the 
model SyH-SNI-SB at t = 3.75 Myr. Density is in cm“^ and ve¬ 
locity is in units of the reference sound speed (c^ 5xio4 = 33 km 
s“^). Here, the velocity signal charaterizes the real velocity direc¬ 
tion (rather than its absolute value). 

iOl)- Thus the corresponding time scale to cross the disk, 
of about 0.4 Myr, is longer than in the previous model, 
during which a slow shock wave propagates into the disk 
transporting mass, momentum and heating to the ISM 
arising its temperature to about 10^ K (with peaks of 
10^ K). 

However, after the jet breaks out into the halo the 
shock wave that propagates along the plane of the disk 
reaches a steady state, separating the low-density, high- 
temperature region where the jet propagates from the 
rest of the disk, within a distance of 80 pc from the center 
of the galaxy. This region is larger than that obtained 
in the previous simulation with a denser and slower jet 
(see Eig. [20]), but this difference does not affect much the 
evolution of the disk and confirms the results obtained 
in the previous models, as we can see in Eigures|2T]and 
|22j where we compared the mass evolution and the mass 


6. DISCUSSION AND CONCLUSIONS 

In this study our main goal was to understand the 
role of the collimated jet emerging from the center of a 
Seyfert-like galaxy in the development of a broader mas¬ 
sive gas outflow in the nuclear region. In contrast to pre¬ 
vious studies (e.g. iCiotti fc Ostrik^l2()Q7l: iSiiacki et ah! 
I2QQ7I: iGaspari et al.ll2QIIbl: lEabianll2QI2f )T we focused our 
attention on a small region of the host galaxy, character¬ 
ized by a volume of 1 kpc^, and performed high resolu¬ 
tion (1.9 pc) hydro dynamical simulations which helped 
us to distinguish the main driving mechanisms of the gas 
evolution. 

It has long been recognized the importance of a SMBH 
jet in balancing cooling flows, inflating the X-ray cavities 
in the hot gas of galaxy clust er cores, or energizing 
the intergalactic med iu m (e .g., iBoehringer et aP 119931: 
McNamara fc NulsenI |2QQ^ lEalceta-Goncalves et al l 
2QlQbl al: IGaspari et ah! l2Qllbl lal). However, all these 
processes act at large scales, from 10 kpc to 1 Mpc. Less 
explored is the importance of the jet feedback on the 
small scales of the host galaxy. 

^ In the other models investigated here we employed solar metal¬ 
licity. 
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Figure 17. Edge-on logarithmic gas density (upper left panel), temperature (upper right panel), pressure (bottom left panel) and velocity 
(absolute value, vertical component, bottom right panel) distribution for the model SyH-SNI-SB-JET large at t = 3.75 Myr. Density is in 
cm“^, temperature in K, pressure in units of 7x10“^^ erg cm“^ (code units) and velocity in units of the reference sound speed (c^ sxio^ 
= 33 km s“^). 


Considering a rotating disk gas distribution initially in 
equilibrium with the galaxy gravitational potential with 
a given column density, we explored the evolution of this 
distribution in the surrounds of the nuclear region over ^ 
10 Myr, taking into account, either separately or all to¬ 
gether, the effects of star formation, SN explosions, and 
the jet. In general lines, we have found that the for¬ 
mation of turbulent and clumpy outflows in the nuclear 
region of Seyfert-like systems is mainly driven by intense 
star formation. The jet feedback at these scales is less 
relevant. 

In the next paragraphs we discuss our main results and 
their implications for the evolution of the nuclear regions 
of Seyfert galaxies, the formation of a multiphase ISM, 
the gas outflow, and the jet feedback on the host galaxy. 

6.1. Gas mass evolution and outflow collimation 

In all models studied, the gas of the disk is carried out 
in an outflow only in presence of star formation activity. 


regardless of the column density of the host galaxy. Nor 
the regular activity in the bulge (mainly due to SNI ex¬ 
plosions), nor the jet, or both together are able to drive 
gas outflow at the small scales, particularly after the first 
Myr, when the jet breaks through the disk and evolves 
along a low-density and high-temperature channel with 
a thickness of few tens of pc. We should stress, however, 
that we considered a radiative cooling optically thin gas 
and neglected the presence of dust and radiation pres¬ 
sure. Though this could be another driving mechanism, 
as remarked in §I the typical gas column densities (be¬ 
tween 10^^ and ^ 10^^ cm^) and the rapid destruction 
rate by SN shock waves make it unlikely the survival of 
substanti al amounts of dust in Seyfert nu clear regions 
(see, e.g. lMcKeel[l989l: I Jones fc Nnthl[2QTTI) . 0 

Though the analytical steady state approach indicates 

® It is interesting to note that in a recent study, I Gan et aU (1201411 
addressed the problem of the AGN feedback in the nuclear regions 
of elliptical galaxies, and argued that Compton scattering can heat 
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Figure 18. Time evolution of the mass of gas within the whole 
system for the models SyH-SNI-SB (top panel) and SyH-SNI-SB- 
JET (bottom panel) for different vertical velocities. Time is in Myr 
and masses are in units of Mq, logarithmic scale. 
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Figure 19. Time evolution of the mass of gas within the whole 
system for the models SyH-SNI-SB (top panel) and SyH-SNI-SB- 
JET (bottom panel) for different gas temperatures. Time is in Myr 
and masses are in units of M©, logarithmic scale. 


that the energy lost by radiative cooling in the denser 
regions of the disk would be high enough to balance the 
energy injected by the stellar feedback, thus preventing 
the formation of a galactic wind (see §3), the numeri¬ 
cal hydro dynamical simulations have demonstrated that 
this steady state lasts only about 1—2 Myr, after which 
the gas begins to flow out of the disk. In fact, at the 
beginning of the simulations the gas in the disk remains 
unperturbed and homogeneous, and its behaviour follows 
the one described by the analytical solutions in §3. How¬ 
ever, with the energy injection by the stellar feedback, 
shock waves propagate into the ISM and the system be¬ 
comes highly non-steady. Shock heated, low density gas 
emerges from the disk and propagates into the stratified 
gas, transferring mass to heights greater than 200—300 
pc. The average mass transfer rate we obtained is be¬ 
tween 3 and 25 Mq yr“^ for a SFR of 1 Mq yr“^ and 
a disk column density between 10^^ and 10^^ cm“^, and 
between 60 and 90 Mq yr“^ for a SFR of 10 Mq yr“^ 
and a disk column density between 10^^ and 10^^ cm“^. 
The gas outflow ends when about ^50—70% of the mass 
in the nuclear region (within a radius of ^ 300 — 400 
pc, is removed. At larger heights, the impact of the SN 
energy on the gas mass evolution is not so important. 
For instance, at a distance of 500 pc above the disk we 
find that the mass loss rate is only a few Mq yr“^ (for a 
SFR of 1 Mq yr“^) and this value is almost independent 
of the disk column density. 

On the other hand, when considering only the propa¬ 
gation of the SMBH jet through the disk there is only a 
small fraction of gas removal from the core of the galaxy 
(within r ^ 40 pc), but no relevant amount is carried 
away by shock waves and/or by gas heating to generate 
a wind. 

We have evaluated also the energy associated to the 
evolved outflow in the models and found that, in the 
cases without the AGN feedback, the kinetic power of 
the outflow increases to a value between (2 to 4)xl0^^ 
erg s“^, that is, about 6.5% to 13% of the stellar feedback 
luminosity. At the same time, the thermal (or internal) 
energy rate of the gas outflow increases to about 10^^ 
erg s“^, that is, ^ 3% of the stellar feedback luminosity. 
In other words, about 15—20% of the stellar luminosity 
goes to the outflow, and the remaining 80% is radiated 
away. In the models where the jet feedback is included, 
we have obtained essentially the same luminosity values 
and for this reason we can conclude that we are observing 
a SF-driven wind. 

The scenario would change completely in the case of 
an AGN outflow emergin g with very large ope ning an¬ 
gles, as assumed, e.g., in IWagner et al.l ()2013F These 
authors (interested in exploring quasar-mode AGN feed¬ 
back) have considered an AGN accretion disk outflow 
impinging into a very cloudy environment with a large 
opening angle of 30 degrees. In this case they find that 
the outflow can substantially affect the galactic ISM. 
Nevertheless, as in the present simulations, their AGN 
wind is injected within a scale of 2 pc. The only dif- 


the gas in the surroundings of the AGN to temperatures as high as 
10® K and for this reason have taken into account the effects of the 
AGN radiative feedback on the system. In their simulations (which 
are performed for two-dimensional ffows that arise from the center 
with a large opening angle), they find that this heating process 
will dominate only in the central regions over radial extensions of 


about 10 to 30 pc. Similarly, in our simulations with AGN feed¬ 
back, although we neglected explicitly the AGN Gompton heating 
effect, we also find that the temperature grows to ~ 10® K within 
a radial extension up to 20-35 pc due to the AGN jet shock heat¬ 
ing. Therefore, the extension of the infiuence of the AGN heating 
feedback found in both studies is essentially the same. 
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Figure 20. Edge-on logarithmic gas density distribution for the models SyH-SNI-SB-JET light (left panels) and SyH-SNI-SB-JET (right 
panels) at the nuclear region between -0.25 and 0.25 kpc (along R and z directions), at t = 0.4 Myr (upper panels) and t = 3 Myr (bottom 
panels). Distances are given in kpc and densities are in cm“^. 


ference is that in our models the jet emerges very col¬ 
limated, as expected, while in IWagner et alJ (j2Q13l ) the 
AGN disk wind is assumed to emerge largely uncolli- 
mated from the nuclear region. This lateral expansion 
may be due to the wind pressure (which in a quasar¬ 
mode AGN may be dominat ed by radia t ion) overcoming 
the external pressure (e.g., iKind I2QQ3I: iCrenshaw et al.l 
l2QQ3l : lArav et al.ll2QQ5f ). However, as stressed, we expect 
that in the nuclear regions of the Seyferts this quasar¬ 
mode feedback and therefore, the development of strong 
AGN disk winds with large opening angles should be in¬ 
operative. 

In summary, the nuclear galactic gas evolution of a 
Seyfert is almost insensitive to the passage of the AGN 
outflow and, as demonstrated above, it needs an intense 
and more widespread source of energy injection in the 
disk (as a high SFR) to increase the gas pressure and 
drive shock waves and turbulence that generate an ex¬ 
panding, multiphase ISM over ^ 500 pc above the disk. 

As discussed in §2, the soft X-ray emission in 
Seyferts is likely dorn inated by photoionized gas (e.g., 
iGnainazzi et al.l 120091) . Its location, mass and mo¬ 
mentum structure resembles a biconical outflow and 


it has been suggested that this might be an ev i dence 
of an AGN-drive n wind (e.g., iLonginotti et al.l 120081: 
lEvans et al.l l201(il ). Nevertheless, a SB-driven outflow 
may also have a biconical shape above the disk (e.g. 
iMelioli et al.ll20I^ . Given the small scales here simu¬ 
lated it is difficult to reproduce the angle of gas escape, 
but the distribution of the kinetic energy that will be 
discussed in §6.3 provides useful information about the 
morphology of the gas outflow. We should also note that 
a new set of numerical simulations at intermediate scales 
(~ 10 kpc) for these flows is in progress which will allow 
to identify more precisely these morphological features 
(Melioli et al. 2015, in prep.). 

6.2. Multiphase ISM 

The SN explosions and the development of shock waves 
into the ISM of the host galaxy determine the formation 
of a multiphase ambient where gas with different densi¬ 
ties and temperatures coexists at a given pressure. In¬ 
deed, almost independently of the presence of a SMBH 
at the center of the galaxy, the stellar feedback gener¬ 
ates naturally regions at high densities which cool fast, 
and regions where the passage of shock fronts leaves low 
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Figure 21. Time evolution of the mass of the gas within the 
whole system (z < 500 pc, upper panel), the thick disk (z < 200 
pc, middle panel) and the central core of the galaxy (r < 40 pc, 
bottom panel) for the models SyH-SNI-SB-JET light (solid-black 
lines) and SyH-SNI-JET light (dotted-blue lines). Time is in Myr 
and mass is in units of M©, logarithmic scale. 



t (Myr) 


Figure 22. Time evolution of the gas mass transfer and loss rate 
of the thick disk (z< 200 pc) for the models SyH-SNI-SB-JET 
light (solid-black lines) and SyH-SNI-JET light (dotted-blue lines). 
Time is in Myr and mass loss rate is in units of Mq yr“^. 


density and very high temperature gas tracks. Thus a 
region characterized by a high SFR may inject enough 
energy to blow-out the gas from the disk and at the same 
time lead to the formation of clouds and clumps which 
can then be continuously generated and steadily carried 
away roughly preserving the total cloud number. 

At the beginning of its evolution, right after launch¬ 
ing, also the jet generates high density regions through 
a transverse shock wave that propagates radially until 
reaching an equilibrium position at about 20—30 pc away 
from the galaxy center (see Figure [20l right panels). This 
shock wave ini tially induces some turbulence (as sng- 
gested by, e.g., iGuillard et al.l[ 20 T^ lAlatalo et al.ll20f5h 
that may trigger either star formation or quenching de¬ 
pending on the shock strength and the ISM conditions 
(see, e.g.. lMelioli et al.l[2006ULeao et al.ll2009[ ). but com¬ 
pared to the dynamical time scales of the system this 
strong interaction lasts too little (about 0.3 Myr) in our 
simulations and it is hard to predict whether or not it 
may have a major impact over the star formation his¬ 
tory of the inner region of the galaxy. 0 

Clouds and clumps are therefore continuously and 
mostly formed in the nuclear region of the galaxy by 
the fragmentation of the shocked gas compressed by the 
supernova shock fronts, and we verify that in our simu¬ 
lations about 30% of the gas mass goes to clouds char¬ 
acterized by high density, between 100 and 1000 cm“^. 
These clouds, that are photoionized by the central source 
(and by the stellar feedback too) and radially dragged out 
of the nuclear region mainly via ram-pressure, could be 
identified with the NLR observed in the Seyfert galaxy. 
On the other hand, we have found that the jet feedback 
alone is unable to produce a multiphase ambient and the 
dense clouds in the nuclear region of the host galaxy. 

At a given distance from the source of high energy 
photons, these clouds may cool very fast to a tempera¬ 
ture of few tens of K, since the cooling time is tcooi ^ 
50/(A(T, 2 ^) 10 - 2471103 ) yr, where A(T, 2 ^) 10-24 is the value 
of the cooling function in units of 10 “^^ erg cm^ s“^ and 
77 io 3 is the gas density of the clouds in units of 10 ^ cm“^. 
From our simulations we note that these clouds are mov¬ 
ing at velocities between -30 and 300 km s“^, relative to 
the galaxy disk, depending on the height above the disk, 
and are characterized by a density contrast between 100 
and 10^, as shown in Figure [23l 

The fact that some clouds have negative velocities im¬ 
plies that they may eventually fall back into the disk in a 
proc ess that may resernble th at of the galactic fountains 
(e.g. lMelioli et al.l[2008Ll200^ . while the clouds with pos¬ 
itive velocities larger than the escape velocity can con¬ 
tinue their journey through the halo. If dust forms in the 
way, it may eventually also originate molecular clouds, as 
recently observed, for exampl e, in the western radio lobe 
of the Seyfert galaxy IC 5063 (|Morganti et al.l[2Q13f ). We 
note however, that in this system they detected clouds 
with higher velocities than those found in our simula¬ 
tions (which are around 1200 km s“^). This may be an 

^ We should note that, in the case (not considered here) of an 
intermittent jet, this would naturally allow for a continuous repe¬ 
tition of events such as those described above, i.e., every time that 
the jet would resume, the digging of a new tunnel into the ISM 
would induce a new shock wave, turbulence and the boosting or 
the quenching of star formation, over the whole evolution of the 
galaxy or as long as the intermittent jet would survive. 
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indication that the observed clouds in this case could be 
produced in the SB-driven wind, as described above, and 
then dragged by the much faster jet flow. Another pos¬ 
sibility is that they result from direct interactions of the 
jet head front with the dense environment. We actually 
see this interaction in the very beggining of the simula¬ 
tions (in the first 10^ yr) when the jet breaks into the 
ISM for the first time and sweeps the interstelar matter 
(see also, e.g., Falceta-Goncalves et ah, in prep.). 

6.3. Velocity distribution 

The velocity distributions depicted in Figures [18] 
clearly show that one of the main differences between 
models with and without a SMBH jet is the gas with 
velocity signatures larger than 1500 km s“^. Depending 
on the column density of the disk of the host galaxy, the 
amount of very high velocity gas when a jet is present 
is between 5 and 20 times larger than in models without 
a jet. However, the fastest moving gas is concentrated 
around the jet and the gas mass with these very high 
velocities corresponds to about 500 — 1000 Mq only, an 
amount too low to be able to modify the overall evolu¬ 
tion of the system, t hough it may influence in the drag of 
high velocity clouds (|Morganti et al.ll^OlSUTombesi et al.l 
I2011) . as remarked in §6.2. 

Therefore, we may conclude that the jet plays an im¬ 
portant role in producing the very high gas velocities 
observed in the expanding nuclear regions of the Seyfert 
galaxies, but at the same time these velocities involve 
such a negligible amount of gas mass compared to the 
total amount in the outflow that it is unable to drive the 
gas exchange between the disk and the halo within the 
kpc scale. 

According to recent observations (see §. 2 and ref¬ 

erences therein), we note that our results basically re¬ 
produce the structures of the gas outflow, that is, an 
extended gas outflow with systemic velocity around the 
nucleus having a broader biconical component perpen¬ 
dicular to the disk and an inner component due to the 
interaction between the jet and the galactic disk mate¬ 
rial. 

Our simulations suggest that in order to achieve this 
structure, three different contributing factors are neces¬ 
sary: i) a diffuse region characterized by high SFR (or 
by a SB) which is responsible for the increase of ther¬ 
mal pressure, the driving of shock waves and the forma¬ 
tion and acceleration of the denser features; ii) a SMBH 
jet, responsible for a little fraction of very high velocity 
clouds and for the gas at very high temperatures (see Fig. 
ttS); and Hi) the presence of high density gas within the 
thin disk (of the order of 10^ cm“^) in order to constrain 
most of the SB driven gas outflow to the inner radii (since 
the SB activity is smaller in the outer radii) and allow 
the formation of the biconical shape as observed. 

Figure 123] that depicts the edge-on logarithmic distri¬ 
bution of the kinetic energy of the gas for models SyH- 
SNI-JET and SyH-SNI-SB-JET large (which differ only 
by the extension of the SB region), illustrates the sce¬ 
nario above. The gas along the jet is clearly at very 
high velocity, but its high collimation makes negligible 
the amount of kinetic energy that is transported from 
the disk to the halo at hz ^ 500 pc. At the same time, 
the presence of a larger SB region causes an almost uni¬ 
form expansion of the gas within 300—500 pc from the 


center of the galaxy. In some regions the flow carries an 
amount of kinetic energy about 1000 times larger than 
that associated to the jet propagation only. This energy 
transfer causes the development of a nearly biconical gas 
outflow characterized by a width of ^ 200 pc and an 
opening angle of ^ 60° at hz = 400 pc. 

We should also notice that several earlier studies have 
ruled out completely SE, or more specifically, SNe as 
the main source of wind feedback based on the high ob- 
serv ed outflow speeds, la r ger than or y ound 600 km/s 
(e.g iMaiolino et al.l [20121: iFabianl r2012[) . However, our 
simulations indicate speeds between 700 and 1500 km/s 
even for models with no AGN feedback (see SyH-SNI-SB 
large case in Eig. nil). As a matter of fact an individual 
SN-driven gas outflow is expected to have velocities of a 
few hundred km/s. However, these velocities are related 
to the colder, denser part the gas that is pushed by the 
external shock shell of the expanding bubble inflated by 
the SN itself. In a more complex system, like the one 
we are investigating here, the stellar feedback generates 
several SN bubbles that interact which each other giving 
rise to a complex multi-phase ambient formed by low ve¬ 
locity, high density gas mixed with low density, very high 
temperature gas moving at much higher speeds. The his¬ 
togram of the gas density distribution versus velocity, in 
Figure [25l illustrates this point for the jetless model SyH- 
SNI-SB large (upper diagram). We clearly see that the 
very small fraction of gas that is accelerated to velocities 
of about 1000 km/s has densities between 10“^ and 10“^ 
cm“^. In other words, the highest outflow speeds are re¬ 
lated to a little fraction of gas at very low density and 
very high temperature. In the same Figure (bottom dia¬ 
gram), we see that the introduction of the jet increases a 
little the gas fraction in the high velocity tail (to values 
larger than 10^ km s“^), but the remainder of the gas 
behaves similarly as the jetless model (see also Figures 
M and nzi bottom-right panel). 

6.4. Conclusions 

In summary our results indicate that: 

• The nuclear galactic gas evolution of a Seyfert is 
almost insensitive to the passage of the jet and it 
needs an intense and more widespread source of en¬ 
ergy injection in the disk to increase the gas pres¬ 
sure and drive shock waves and turbulence; 

• The SN explosions and the development of shock 
waves into the ISM of the host galaxy determine 
the formation of a multiphase ambient, and there¬ 
fore a SB and/or a region with high SFR may in¬ 
ject enough energy to blow-out the gas from the 
disk and at the same time lead to the formation 
of clouds and clumps which can be continuously 
generated and steadily carried away. Therefore, 
we may roughly expect a preservation of the total 
cloud number; 

• The fact that some clouds have negative velocities 
implies that they may eventually fall back into the 
disk, while the clouds with positive velocities larger 
than the escape velocity can continue their journey 
through the halo. If dust forms in the way, it may 
eventually also originate molecular clouds; 
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Figure 23. Face-on logarithmic gas density (top-left panel), temperature (top-right panel), pressure (bottom-left panel) and velocity 
(bottom-right panel) distribution at t = 6.5 Myr, for model SyH-SNI-SB-JET large at a height h = 430 pc. Distances are given in kpc, 
densities, temperatures and pressures are in cgs and velocities are in units of the reference sound speed, computed at T=5xl0^ K, c^ 

= 33 km s“^. 


• The velocity distribution reproduced by our simu¬ 
lations clearly show that one of the main differences 
between models with and without a SMBH jet is 
the presence of gas with velocity signatures larger 
than 1500 km s“^. Depending on the column den¬ 
sity of the disk of the host galaxy, the amount of 
very high velocity gas when a jet is present is be¬ 
tween 5 and 20 times larger than in models without 
a jet; 

• The jet plays an important role in producing the 
very high gas velocities observed in the expand¬ 
ing nuclear regions of the Seyfert galaxies, but at 
the same time these velocities involve such a neg¬ 
ligible amount of gas mass (compared to the total 
amount in the outflow) that it is unable to drive 
the gas evolution within the kpc scale. Neverthe¬ 
less, they may be responsible for recent observed 
high velocity feature s in Seyfert and ULIR Gs (e.g. 
iMorganti et al.|[2QTl: iTombesi et al.|[2QT5l) : 


• Our results basically reproduce the structures of 
the gas outflow, that is, an extended gas outflow 
with systemic velocity around the nucleus having a 
broader biconical component perpendicular to the 
disk and an inner component due to the interaction 
between the jet and the galactic disk material. 

Finally, we should remark that we have neglected 
the effects of magnetic fields in the present study. 
Although they are known to be particularly im¬ 
portant in the formation and propagation of the 
jet outflow, we have neglected their effects here for 
simplicity. Nevertheless, we may expect that its 
presence would not affect much the overall conclu¬ 
sions of this work. In fact, the inclusion of mag¬ 
netic fields in the shaping of the jet would help 
on one side to confine the beam even more in the 
central narrow channel, therefore further constrain¬ 
ing its interaction with the surroundings. On the 
other side, the enhanced collimation and magneto- 
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Figure 24. Edge-on logarithmic kinetic energy gas distribution 
for models SyH-SNI-SB-JET large (left panel) and SyH-SNI-JET 
(right panel) at t = 4.8 Myr. Distances are given in kpc and energy 
is in units of 5x10"^^ erg (code units) 


centrifugal acceleration implied by the presence of 
magnetic fields in the jet would help to push cap¬ 
tured clumps from the underlying star-formation 
driven wind to much higher speeds as those sug¬ 
gested by the recent observations mentioned above. 
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Figure 25. Two dimensional histogram of the vertical velocity 
versus density 

calculated considering every cell within the whole simulated system 
for the models SyH-SNI-SB large (top panel) and SyH-SNI- 
SB-JET large (bottom panel) at a time t=5 Myr. Velocities 
are in units of km s“^ (logarithmic scale) and densities are in units 
of cm“^ (logarithmic scale). The colour bar indicates the cell 
number normalized to their total number. 

group of lAG-USP. 


APPENDIX 

A. SEYFERT GALAXY MODELS WITH LOW GOLUMN DENSITY (SYL) 

We have also run numerical models considering the presence (or not) of SNI, SB and JET in a disk galaxy charac¬ 
terized by lower values of the column density (SyL). The results are consistent with those presented above, and the 
gas evolution in the disk is completely driven by the SB, as can be seen in Eig. [26l where the gas density evolution 
for each model (SyL-SNI, SyL-SNI-SB, SyL-SNI-SB-JET and SyL-SNI-JET) at t=1.5 Myr is shown. 

Examining Eigure[23we note that when a SB region is active, regardless of the presence of a SMBH jet, the mean 
























23 



Figure 26. Edge-on logarithmic gas density distribution at t = 1.5 Myr, for the models SyL-SNI (top-left panel), SyL-SNI-SB (top-right 
panel), SyL-SNI-SB-JET (bottom-left panel) and SyL-SNI-JET (bottom-right panel). Distances are given in kpc and densities are in cm“^ 


mass transfer rate from the the disk is about 3 Mq yr“^ over a time of ^2 Myr, resulting in a net total gas mass lost 
by the system of about 2 xlO^ Mq. Also the velocity distribution, shown in bottom-right panel of Fig. [271 confirms 
the great similarity between the two models with and without the SMBH jet. Both models lost about 40% of their 
mass, and have about 95% of the remaining mass characterized by the same velocity profile. 
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